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ABSTRACT 

Motivated  by  the  recent  interest  in  ocean  energetics,  the  widespread  use  of  horizontal  eddy  viscosity  in 
models,  and  the  promise  of  high  horizontal  resolution  data  from  the  planned  wide-swath  satellite  altimeter, 
this  paper  explores  the  impacts  of  horizontal  eddy  viscosity  and  horizontal  grid  resolution  on  geostrophic 
turbulence,  with  a  particular  focus  on  spectral  kinetic  energy  fluxes  n(/Q  computed  in  the  isotropic  wave- 
number  ( K)  domain.  The  paper  utilizes  idealized  two-layer  quasi  geostrophic  (OG)  models,  realistic  high- 
resolution  ocean  general  circulation  models,  and  present-generation  gridded  satellite  altimeter  data.  Adding 
horizontal  eddy  viscosity  to  the  QG  model  results  in  a  forward  cascade  at  smaller  scales,  in  apparent 
agreement  with  results  from  present-generation  altimetry.  Eddy  viscosity  is  taken  to  roughly  represent 
coupling  of  mesoscale  eddies  to  internal  waves  or  to  submesoscale  eddies.  Filtering  the  output  of  either  the 
QG  or  realistic  models  before  computing  n(/Q  also  greatly  increases  the  forward  cascade.  Such  filtering 
mimics  the  smoothing  inherent  in  the  construction  of  present-generation  gridded  altimeter  data.  It  is 
therefore  difficult  to  say  whether  the  forward  cascades  seen  in  present-generation  altimeter  data  are  due  to 
real  physics  (represented  here  by  eddy  viscosity)  or  to  insufficient  horizontal  resolution.  The  inverse  cascade 
at  larger  scales  remains  in  the  models  even  after  filtering,  suggesting  that  its  existence  in  the  models  and  in 
altimeter  data  is  robust.  However,  the  magnitude  of  the  inverse  cascade  is  affected  by  filtering,  suggesting  that 
the  wide-swath  altimeter  will  allow  a  more  accurate  determination  of  the  inverse  cascade  at  larger  scales  as 
well  as  providing  important  constraints  on  smaller-scale  dynamics. 
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1.  Introduction 

This  paper  examines  the  impacts  of  horizontal  eddy 
v  iscosity  and  of  horizontal  grid  resolution  on  geostrophic 
turbulence.  Particular  emphasis  is  placed  on  spectral 
fluxes  U(K)  of  surface  ocean  geostrophic  kinetic  energy 
computed  in  the  isotropic  wavenumber  (K)  domain  from 
an  idealized  two-layer  quasigeostrophic  (QG)  turbulence 
model,  from  a  high-resolution  ocean  model  run  in 
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a  realistic  near-global  domain  (hereafter,  the  “realistic” 
model),  and  from  present-generation  gridded  satellite 
altimeter  data.  We  are  motivated  by  considerations  of 
ocean  energetics  and  by  the  widespread  use  of  horizontal 
eddy  viscosities  in  models.  We  are  especially  interested  in 
whether  forward  kinetic  energy  cascades  may  represent  a 
significant  sink  of  energy  in  the  ocean  and  may  be  trig¬ 
gered  by  processes  that  can  be  roughly  characterized  by 
a  horizontal  eddy  viscosity.  The  likelihood  of  greatly  in¬ 
creased  horizontal  resolution  in  planned  future-generation 
wide-swath  satellite  altimeter  data  (Fu  and  Ferrari  2008; 
Fu  et  al.  2012)  motivates  us  to  examine  here  the  im¬ 
portance  of  horizontal  grid  resolution  in  the  estimation 
of  spectral  fluxes,  particularly  at  smaller  scales.  An  ad¬ 
ditional  goal  of  the  paper  is  to  directly  compare  the  in¬ 
verse  cascades  at  larger  scales  diagnosed  in  spectral 
fluxes  n(/Q  in  realistic  models  with  the  inverse  cascades 
seen  in  present-generation  gridded  altimeter  data,  as  a 
test  of  the  robustness  of  both.  The  model-data  compari¬ 
son  will  be  conducted  over  six  oceanic  regions,  which  vary 
widely  in  their  eddy  kinetic  energy  levels. 

Recently  the  oceanographic  community  has  focused 
considerable  attention  on  quantifying  the  mechanisms 
and  spatial  distribution  of  energy  input,  energy  dissipa¬ 
tion,  and  eddy  diffusivity  in  the  ocean.  The  recent  interest 
in  ocean  energetics  has  arisen  in  part  because  ocean 
mixing  is  thought  to  exert  a  strong  control  on  abyssal 
stratification  and  circulation  (e.g.,  Munk  and  Wunsch 
1998).  One  important  energy  source  is  the  ~1-TW  wind 
power  input  into  geostrophic  flows  (e.g.,  Wunsch  1998; 
Scott  1999;  Scott  and  Xu  2009).  Several  recent  papers 
have  explored  potential  dissipation  mechanisms  for  this 
energy.  Sen  et  al.  (2008),  Arbic  et  al.  (2009),  Wright  et  al. 
(2013),  and  references  therein  argue  that  quadratic  bot¬ 
tom  boundary  layer  drag  is  an  important  energy  sink  for 
low-frequency  flows.  Marshall  and  Naveira-Garabato 
(2008),  Nikurashin  and  Ferrari  (2011),  and  Scott  et  al. 
(2011)  argue  that  internal  lee  waves  generated  by  geo¬ 
strophic  flows  over  rough  topography  are  a  significant 
energy  sink.  An  additional  mechanism  for  dissipation 
over  rough  topography  has  recently  been  proposed  by 
Dewar  and  Hogg  (2010). 

Here  we  consider  another  potential  mechanism,  namely, 
a  forward  cascade  of  kinetic  energy  to  small  scales,  partic¬ 
ularly  in  the  upper  1000  or  so  meters  of  the  ocean  (i.e.,  the 
thermocline).  In  the  models  used  in  this  paper,  the  forward 
kinetic  energy  cascade  will  be  triggered  and  dissipated 
by  horizontal  eddy  viscosity  v  (hereafter,  “viscosity”).1 
Viscosity  has  long  been  employed  in  numerical  ocean 


1  The  important  subject  of  vertical  eddy  viscosity  will  not  be 
addressed  in  this  paper. 


models  (e.g.,  Muller  1979,  among  many).  Viscosity  is  used 
to  absorb  subgrid-scale  structure  in  models  and  is  generally 
kept  as  small  as  the  model  resolution  allows  (e.g.,  Wallcraft 
et  al.  2005).  The  effects  of  viscosity  on  ocean  climate  models 
have  been  investigated  by  Jochum  et  al.  (2008).  To  our 
knowledge,  this  paper  presents  the  first  study  of  the  impacts 
of  viscosity  on  the  statistical  properties  of  eddies,  including 
spectral  kinetic  energy  fluxes  I1(X),  in  idealized  QG  tur¬ 
bulence  models.  Although  the  focus  of  the  QG  sensitivity 
study  is  on  the  spectral  kinetic  energy  fluxes  n(X),  wc  will 
also  briefly  discuss  the  sensitivity  of  the  energy  levels,  hor¬ 
izontal  scales,  and  vertical  structure  of  eddies  in  the  QG 
turbulence  model  to  viscosity.  We  will  examine  the  impact 
of  viscosity  on  the  energy  budget  of  the  QG  model.  The 
sensitivity  study  to  viscosity  presented  here  serves  as 
a  complement  to  our  earlier  studies  of  the  sensitivity  of  QG 
turbulence  to  linear  bottom  drag  (Arbic  and  Flierl  2004; 
Arbic  et  al.  2007;  see  also  Thompson  and  Young  2006, 
2007).  See  also  Arbic  and  Scott  (2008)  for  an  examination 
of  the  sensitivity  of  QG  turbulence  to  quadratic  bottom 
drag.  In  addition,  we  will  examine  the  impact  of  viscosity  on 
the  spectral  kinetic  energy  fluxes  IT(X)  computed  from 
realistic  models.  We  emphasize  again  that,  for  the  sake  of 
simplicity,  we  focus  here  on  the  spectral  fluxes  of  geo¬ 
strophic  kinetic  energy.  Any  ageostrophic  motions  that  may 
be  present  in  the  realistic  model  will  be  removed,  at  least 
roughly,  with  a  low-pass  filter  in  the  time  domain. 

The  spectral  fluxes  II(X)  reveal  the  direction  of  the 
nonlinear  energy  cascade  and  were  computed  from 
present-generation  gridded  satellite  altimeter  data  by 
Scott  and  Wang  (2005),  following  similar  pioneering 
calculations  performed  on  atmospheric  data  by  Boer 
and  Shepherd  (1983).  Scott  and  Arbic  (2007)  examined 
spectral  kinetic  energy  fluxes  IT(/C)  in  the  baroclinic 
mode,  barotropic  mode,  and  upper  layer  of  a  two-layer 
QG  model  and  used  the  model  results  to  interpret  the 
inverse  cascade  of  kinetic  energy  seen  in  altimetry  data 
by  Scott  and  Wang  (2005)  at  larger  scales.  Scott  and 
Arbic  (2007)  did  not  explore  the  forward  cascade  seen 
in  altimetry  data  at  smaller  scales.  Here  we  focus  much 
more  attention  on  this  forward  cascade.  One  hypoth¬ 
esis  explored  in  this  paper  is  that  the  forward  cascade 
seen  in  present-generation  altimeter  data  is  consistent 
with  the  action  of  processes  that  can  be  roughly  char¬ 
acterized  by  a  horizontal  eddy  viscosity. 

We  will  not  discuss  in  great  detail  here  the  mecha¬ 
nisms  that  may  lie  behind  an  oceanic  forward  kinetic 
energy  cascade.  However,  a  few  brief  remarks  are  in 
order.  Capet  et  al.  (2008)  demonstrated  that  fronto- 
genesis  in  regions  of  active  submesoscale  eddies  yields 
a  forward  kinetic  energy  cascade  that  drains  energy  from 
the  mesoseale.  The  forward  cascade  in  the  Capet  et  al. 
(2008)  results  is  driven  fundamentally  by  ageostrophic 
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velocities.  In  support  of  the  frontogenesis  mechanism, 
D’Asaro  et  al.  (2011)  presented  observational  evidence 
of  enhanced  energy  dissipation  in  ocean  fronts.  Another 
postulated  mechanism  for  a  forward  cascade  is  inter¬ 
actions  between  mcsoscale  eddies  and  the  internal  wave 
field.  Building  upon  studies  of  wave-mean  flow  inter¬ 
actions  (Muller  and  Olbers  1975;  Muller  1976;  Blihler  and 
McIntyre  2005),  Polzin  (2008, 2010)  argued  that  mesoscale 
eddy-internal  wave  interactions  can  be  parameterized 
by  a  viscosity.  Based  on  an  updated  analysis  of  Polygon 
Mid-Ocean  Dynamics  Experiment  (POLYMODE)  data, 
Polzin  (2008,  2010)  finds  comparable,  but  somewhat 
smaller,  estimates  of  viscosities  than  Brown  and  Owens 
(1981)  obtained  from  the  same  dataset. 

We  will  also  not  contribute  here  to  recent  discussions 
(e.g.,  LaCasce  2012  and  references  therein)  about  whether 
the  altimeter  signal  consists  primarily  of  the  first  baroclinic 
mode  of  an  “interior  QG”  solution  (Wunsch  1997),  or  of 
a  “surface  quasi-geostrophic”  (SQG)  solution  (Lapeyre 
2009).  The  idealized  two-layer  QG  model  used  in  our 
study  does  not  include  any  SQG  dynamics,  and  the  re¬ 
alistic  model  employed  here  does  not  have  enough 
vertical  resolution  to  properly  include  SQG  effects. 
Our  choice  of  models  is  made  based  on  what  is  easily 
available  to  us,  and  not  because  we  believe  SQG  dy¬ 
namics  to  be  unimportant. 

Motivated  by  the  relatively  coarse  resolution  of 
present-day  gridded  satellite  altimeter  data,  this  paper 
addresses  the  impact  of  horizontal  grid  resolution  on  the 
estimation  of  spectral  kinetic  energy  fluxes.  Along-track 
altimeter  data,  which  is  one-dimensional,  cannot  be  used 
for  spectral  flux  computations,  which  require  a  two- 
dimensional  field.  Instead,  observationally  based  com¬ 
putations  of  spectral  fluxes  must  rely  on  two-dimensional 
gridded  altimeter  products.  The  Archiving,  Validation, 
and  Interpretation  of  Satellite  Oceanographic  (AVISO) 
gridded  altimeter  data  (Le  Traon  et  al.  1998;  Ducet  et  al. 
2000)  used  by  Scott  and  Wang  (2005)  to  compute  spectral 
kinetic  energy  fluxes  are  constructed  by  interpolating 
along-track  data  from  two  satellite  altimeters,  in  space 
and  time,  onto  a  V30  Mercator  grid,  with  7-day  spacing  in 
time.  The  Ocean  Topography  Experiment  (TOPEX)/ 
Poseidon  and  Jason  altimeters  utilized  in  this  procedure 
have  high  along-track  resolution  (about  6  km),  but  a 
large  distance  between  adjacent  tracks  and  a  long  sam¬ 
pling  period  (about  300  km  and  9.9  days,  respectively; 
Chelton  et  al.  2001).  Because  of  the  large  cross-track 
distance  and  long  sampling  period,  the  feature  resolu¬ 
tion  of  gridded  AVISO  data  is  much  coarser  than  6  km 
(Chelton  et  al.  2011).  The  interpolation  procedure  used 
to  construct  gridded  AVISO  data  acts  to  smooth  the 
raw  along-track  data,  in  both  space  and  time.  Motivated 
by  this  smoothing,  we  examine  here  spectral  kinetic 


energy  fluxes  II(X)  estimated  from  filtered  versus  unfil¬ 
tered  output  of  our  models. 

Capet  et  al.  (2008)  noted  that  the  zero  crossings 
(wavenumbers  at  which  the  spectral  kinetic  energy  fluxes 
U(K)  cross  zero,  i.e.,  change  sign  from  negative  to  posi¬ 
tive)  were  rather  different  in  their  simulations  compared 
to  the  Scott  and  Wang  (2005)  altimetrically-derived  fluxes. 
Capet  et  al.  (2008)  noted  that  “it  is  unclear  if  the  forward 
cascade  present  in  these  [altimeter]  observations  occurs 
for  the  same  reasons  than  in  our  [the  Capet  et  al.  (2008)] 
simulations.”  Our  study  of  the  effects  of  smoothing  on 
n(X)  estimates  will  shed  some  light  on  whether  these 
discrepancies  may  be  due  at  least  in  part  to  sampling 
artifacts.  In  effect,  then,  our  alternative  hypothesis  for 
the  forward  cascades  seen  in  present-generation  al¬ 
timeter  data  is  that  they  are  due  to  smoothing  in  the 
AVISO  data.  In  this  exercise  we  are  also  m6tivated  by 
Boer  and  Shepherd  (1983),  who  discussed  the  possi¬ 
bility  that  the  forward  cascades  they  saw  at  small  scales 
were  artifacts  of  the  resolution  limitations  in  the  at¬ 
mospheric  reanalyses  they  used. 

The  planned  wide-swath  altimeter  will  have  repeat 
periods  comparable  to  those  of  TOPEX/Jason;  22  days 
has  been  discussed  as  a  target,  for  example  (Fu  et  al. 
2012).  As  we  will  show,  temporal  smoothing  can  distort 
the  spectral  fluxes  n(X)  of  geostrophic  kinetic  energy 
computed  from  the  sea  surface  height  (SSH)  field.  How¬ 
ever,  because  snapshots  of  wide-swath  data  will  be  two- 
dimensional  (and  at  high  horizontal  resolution),  it  will  be 
possible  to  compute  spectral  fluxes  II(X)  directly  from  the 
snapshots.  Since  the  spatial  and  temporal  smoothing 
needed  to  construct  current-generation  gridded  altimeter 
datasets  will  be  much  reduced,  spectral  fluxes  n(AT) 
computed  from  wide-swath  altimetry  should  be  much 
more  reliable.  Prompted  by  a  reviewer,  we  emphasize 
here  that  altimetry,  even  high-resolution  wide-swath  al¬ 
timetry,  can  only  be  used  to  compute  spectral  fluxes  II(X) 
arising  from  geostrophic  velocities.  Spectral  fluxes  arising 
from  ageostrophic  effects  can  be  computed  from  models 
(e.g.,  Capet  et  al.  2008),  whose  validity  at  smaller  scales 
can  be  checked  by  comparison  of  the  model  SSH  fields 
with  the  wide-swath  altimeter  data. 

We  will  compare  spectral  fluxes  II(/Q  computed  from 
fully  resolved  outputs  of  the  idealized  QG  model  and  of 
the  realistic  model,  with  spectral  fluxes  computed  from 
spatially-  and/or  temporally  filtered  versions  of  these 
outputs.  The  fluxes  computed  from  filtered  versions  of 
the  models  are  taken  to  roughly  represent  the  state-of- 
the-art  in  regards  to  what  can  be  done  with  present-day, 
somewhat  coarse-resolution  gridded  altimeter  products. 
The  fluxes  computed  from  unfiltered  model  output  are 
taken  to  roughly  represent  what  could  be  done  with  the 
proposed  wide-swath  satellite  altimeter.  We  will  focus 
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on  the  impact  of  horizontal  resolution  on  the  estimation 
of  the  forward  cascade  at  small  scales,  but  will  examine 
resolution  effects  on  the  estimated  inverse  cascade  at 
larger  scales  as  well. 

The  paper  is  organized  as  follows.  Section  2  orients 
the  reader  by  presenting  a  typical  spectral  flux  profile 
computed  from  present-generation  gridded  altimeter  data. 
Section  3  presents  a  description  of  the  QG  turbulence 
model,  a  discussion  of  the  impacts  of  viscosity  on  the  QG 
model  behaviors,  and  a  discussion  of  the  impact  of  both 
spatial  and  temporal  filtering  on  the  spectral  fluxes  n(/Q 
computed  from  the  QG  model.  Section  4  presents  spectral 
fluxes  U(K)  computed  from  unfiltered  and  filtered  output 
of  the  realistic  model,  alongside  fluxes  computed  from 
present-generation  gridded  satellite  altimeter  (AVISO) 
data.  A  summary  and  discussion  are  presented  in  section  5. 

2,  Example  spectral  flux  curve  from  gridded 
altimeter  data 

To  set  the  stage  for  discussions  to  follow,  in  Fig.  1  we 
display  the  spectral  fluxes  II(/Q  computed  over  the 
Kuroshio  region  from  AVISO  gridded  altimeter  data. 
Details  of  the  AVISO  II(/Q  calculations  will  be  dis¬ 
cussed  later,  in  section  4.  Figure  1  illustrates  the  nearly 
universal  shape  seen  in  the  spectral  kinetic  energy  fluxes 
n(K)  computed  throughout  the  Pacific  Ocean  by  Scott 
and  Wang  (2005).  A  positive  lobe  at  smaller  scales 


FIG.  1 .  Spectral  fluxes  fI(/Q  of  surface  ocean  geostrophic  kinetic 
energy  vs  isotropic  wavenumber  K  in  Kuroshio  region,  computed 
from  AVISO  gridded  satellite  altimeter  data  as  described  later  in 
section  4.  Dashed  vertical  line  denotes  deformation  wavenumber 
1  ILd  for  this  region.  The  data  source  for  the  Ld  value  is  also  de¬ 
scribed  later  in  section  4. 

(larger  wavenumbers)  indicates  a  forward,  or  direct, 
cascade  of  energy  toward  yet  smaller  scales.  A  larger 
negative  lobe  at  larger  scales  indicates  an  inverse  cas¬ 
cade  of  kinetic  energy  toward  yet  larger  scales.  The 
relative  size  of  the  forward  and  inverse  cascades  can  be 
measured  by  the  quantity  Ftol,  which  we  define  by 


maximum  positive  value  of  U(K)  at  smaller  scales 
minimum  negative  value  of  Yl(K)  at  larger  scales 


a) 


For  the  II(/Q  curve  in  Fig.  1,  the  Ftol  value  is  0.19,  which 
lies  in  the  middle  of  the  Ftol  values  computed  later  in 
section  4  over  six  different  oceanic  regions.  The  de¬ 
formation  wavenumber  1  !Ld,  where  Ld  is  the  first  baro- 
ciinic  mode  deformation  radius,  is  shown  as  a  dashed 
vertical  line  in  Fig.  1.  The  data  source  for  the  oceanic  Ld 
values  used  in  Fig.  1  and  elsewhere  in  the  paper  will  also  be 
described  in  section  4.  For  the  U(K)  curve  shown  in  Fig.  1, 
the  deformation  wavenumber  lies  at  the  zero  crossing — 
the  transition  between  the  negative  and  positive  lobes. 
Later  we  will  see  that  the  transition  wavenumber  in  the 
models  we  examine  is  higher  than  1  !Ld  unless  the  models 
employ  relatively  large  horizontal  eddy  viscosity. 

3.  The  idealized  quasigeostrophic  model 

a.  Description  of  idealized  model 

We  use  the  same  idealized  two-layer  QG  turbulence 
model  as  in  many  of  our  previous  papers  (e.g.,  Arbic  and 


Flierl  2004;  see  also  Flierl  1978).  The  model  is  forced 
by  an  imposed  mean  flow,  which  is  taken  to  roughly 
represent  the  time-mean  flows  of  a  subtropical  gyre.  The 
mean  flow  is  horizontally  homogeneous  within  each 
layer.  The  model  is  doubly  periodic,  is  run  here  on  an 
/  plane  with  a  flat  bottom  for  simplicity,  and  has  256 
gridpoints  on  a  side  (thus  2562  gridpoints  in  all),  where 
each  side  has  dimensional  length  20nLd  except  where 
noted.  The  mean  flow  is  vertically  sheared  between  the 
layers  and  is  therefore  baroclinically  unstable.  The  model 
is  initialized  with  a  randomly  generated  initial  condition 
and  achieves  statistical  equilibrium  when  the  energy  ex¬ 
tracted  by  eddies  from  the  baroclinically  unstable  mean 
flow  is  balanced  by  energy  dissipation.  For  this  paper,  in 
contrast  to  our  earlier  papers  on  QG  turbulence,  we  in¬ 
clude  a  viscosity  in  most  of  our  simulations. 

The  quantity  in  our  two-layer  QG  model  that  most 
closely  corresponds  to  altimeter  measurements  of  SSH 
is  the  upper-layer  streamfunction.  Therefore,  we  focus 
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the  majority  of  our  attention  in  this  paper  on  a  layer 
interpretation  of  the  two-layer  QG  model.  The  gov¬ 
erning  equations  of  the  model  in  layer  form  are 

^L  +  ri  +J(ipvqx)  =  wf  +  uV4ipv  and  (2) 

+  r2 + J(*l>2,  q2)  =  -k2vV2  +  wf + ’  (3) 

di 


where  t  is  time,  potential  vorticity  q  is  the  prognostic 
variable,  and  the  subscripts  1  and  2  denote  the  upper 
and  lower  layers,  respectively.  The  zonal  (east-west) 
and  meridional  (north-south)  spatial  coordinates  are 
respectively  x  and  y,  and  J(A,  B)  =  dAldxdB/dy  —  dA/ 
dydBldx .  Fluctuations  from  the  time  mean  are  written 
without  overbars,  while  imposed  time-mean  quantities 
are  designated  with  overbars.  The  forcing  terms  are 


r  Sql  dqld»l/l 
r«  =  u*  — -  h — - — L 
1  1  ax 


.  1  *  and 


dy  dx  9 


(4) 


2  2  dx  dy  dx  4 


(5) 


The  zonal  velocity  is  w,  the  meridional  velocity  is  v ,  and 
the  imposed  time-mean  flow  is  taken  to  be  zonal.  A 
wavenumber  filter  used  for  subgrid  scale  dissipation  is 
denoted  by  wf  As  shown  by  LaCasce  (1996),  the 
wavenumber  filter  does  not  distort  modons  (Stem  1975; 
Larichev  and  Reznik  1976a, b),  which  are  exact  solu¬ 
tions  of  the  inviscid  nonlinear  equations,  nearly  as  much 
as  horizontal  eddy  viscosities  or  hyperviscosities  do. 
This  suggests  that  the  wavenumber  filter  is  a  very 
“clean”  choice  for  subgrid-scale  dissipation  in  a  model. 
As  discussed  in  Arbic  and  Flierl  (2004),  the  wave- 
number  filter  has  little  effect  (less  than  3%)  on  the 
energy  budgets  of  the  two-layer  forced-dissipated  QG 
turbulence  simulations  in  that  paper.  This  further 
demonstrates  that  the  wavenumber  filter  is  a  “clean” 
choice  for  subgrid-scale  dissipation.  Horizontal  eddy 
viscosity  is  denoted  by  v.  The  perturbation  stream- 
functions  </q  and  if/ 2  satisfy 


<7i  =  VV,  + 


(*A2  ~  i) 

(1  +8)1% 


q2  =  + 


~  t/>2) 
(1  +S)L2d' 


(6) 


where  6  =  H\IH2,  the  ratio  of  upper  to  lower  layer 
depths,  and  the  first  baroclinic  mode  deformation  radius 
Ld  is  defined  as  in  Flierl  (1978).  The  imposed  mean  PV 
gradients  are 

_  s(zq  -  gp 

dy  (1  +  8)1%  dy  (1  +  8)L2' 


The  bottom  Ekman  drag  coefficient  R2  is  determined  by 
the  bottom  boundary  layer  thickness  dekman  =  \/A//o  via 


*2  = 


■^O^Ekman 

2  H2  ’ 


(8) 


where  /0  is  the  Coriolis  parameter  and  A  is  the  vertical 
eddy  viscosity  (cf.  Vallis  2006). 

The  barotropic  (BT)  and  baroclinic  (BC)  stream- 
functions  are  defined  as  in  Flierl  (1978): 


_S*,+*2  _  -  tp2) 

V'bt  fTs-’  V'bc  r+S 


(9) 


The  barotropic  kinetic  energy  (KEBT)  and  baroclinic 
kinetic  energy  (KEbc)  are  given  by 

KEbt=  IJ^btI  2dxdy, 

KEBC=  jj^Bcl2^^-  (1Q) 

The  depth-averaged  energy  equation  in  physical  space 
is  obtained  from  multiplication  of  (2)  by  —6^/(1  4-6), 
multiplication  of  (3)  by  ~if/2l(l  +  6),  integration  over  the 
domain,  and  addition  of  the  results,  yielding: 


l  ffi  [aw2  W2  sQAi  ~ 4>2f] 

arjj2[l+s  1+8  (i+S)2L2J 


dx  dy 


_  3(“l  -  “2) 
(1  +S)2L2 


^(Vip2)2dxdy 


—  V 


dxdy+  wf , 


(id 


where  the  first  term  on  the  right-hand  side  represents 
energy  production,  the  second  term  represents  energy 
dissipation  by  bottom  friction,  the  third  represents  dis¬ 
sipation  by  eddy  viscosity,  and  here  wf  represents  energy 
dissipation  by  the  wavenumber  filter. 

Three  nondimensional  parameters  control  the  be¬ 
havior  of  the  two-layer  QG  model.  Following  Flierl 
(1978)  and  Fu  and  Flierl  (1980),  the  stratification  pa¬ 
rameter  6  is  set  to  a  representative  midlatitude  oceanic 
value  of  0.2  in  all  of  the  runs  here.  The  nondimensional 
linear  bottom  friction  strength  is  =  R2Ldl{u\  —  u2). 
In  simulations  performed  with  viscosity,  the  non- 
dimensional  viscosity  is  vi[Ld(ui  -  u2)].  For  the  sake  of 
convenience  we  will  use  v  to  denote  both  dimensional 
and  nondimensional  viscosity;  the  reader  can  discern 
which  of  the  two  is  meant  by  context. 
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(e)  v=0.1 

20 
0 

-20 

(9)  v-1 

20 
0 

-20 


(f)qi(v=0.1 

■  I  20 

Mo 

1-20 


(h)  q^,  v=1 


Fig.  2.  Snapshots  of  upper-layer  streamfunction  ^  and  upper-layer  potential  vorticity  qx  for 
two-layer  QG  simulations  with  nondimensional  bottom  friction  FL  =  0.4  and  nondimensional 
viscosity  values  v  of  0, 0.01, 0.1 ,  and  1,  respectively.  The  domain  size  is  207 iLd  on  a  side,  where 
Ld  is  the  first  baroclinic  mode  deformation  radius. 
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Figure  2  displays  snapshots  of  the  upper-layer  stream- 
function  if/ 1  and  upper  layer  potential  vortidty  qx  in 
simulations  with  FL  -  0.4  and  nondimensional  v  values  of 
0, 0.01, 0.1,  and  1.  The  v  =  0,  FL  =  0.4  experiment  is  the 
“nominal”  experiment  of  Arbic  and  Flierl  (2004),  that  is, 
the  experiment  having  amplitudes,  horizontal  scales,  and 
vertical  structure  of  eddy  kinetic  energy  that  agree  rea¬ 
sonably  well  with  observations  in  midlatitudes.  The  v  -  0, 
Fl  =  0.4  experiment  consists  of  densely  packed  vorti¬ 
ces  (e.g.,  McWilliams  1984),  in  rough  similarity  to 
a  typical  snapshot  of  altimetrically  measured  SSH  in 
the  midocean  (see,  for  instance,  the  model-altimeter 
comparison  in  Fig.  10  of  Arbic  and  Flierl  2004).  The 
vortices  shown  in  Fig.  2  are  not  as  dominant  as  the 
vortices  in  experiments  with  weaker  bottom  friction 
(see  for  instance  Fig.  8  of  Arbic  and  Flierl  2004).  As  is 
to  be  expected.  Fig.  2  demonstrates  that  increasing 
horizontal  eddy  viscosity  yields  increasingly  smooth 
streamfunction  and  potential  vorticity  fields. 


b.  Spectral  fluxes  in  idealized  model 

As  in  Scott  and  Arbic  (2007),  which  in  turn  followed 
Frisch  (1995),  we  develop  the  spectral  fluxes  II(/Q  by 
letting 

k  t 

+2(x,y,t)  =  I  (12) 

k  t 

where  ^  denotes  a  Fourier  transform  of  ip  in  wave- 
number  space,  and  k  and  /  are  zonal  and  meridional 
wavenumbers,  respectively.  We  neglect  the  subgrid- 
scale  dissipation  w/for  now  since  its  impact  on  spectral 
fluxes  is  difficult  to  compute  directly.  We  will  compute 
its  effects  as  a  residual  of  the  other  terms  in  the  spectral 
flux  equations;  later  we  will  see  that  this  residual  is  small. 
The  Fourier  transformed  upper-  and  lower-layer  gov¬ 
erning  equations  are  then 


1 


I  K2  + 

n 


1 


1 


-s‘k2+ 


+  r2(k,l,t)  =  R2K24>2(k,l,t)  +  vK4<p2(k,l,t), 


(13) 


(14) 


respectively,  where  K2  =  k2  +  l  and  we  have  dropped 
the  summation  signs  for  simplicity. 

We  now  multiply  the  Fourier  transformed  equa¬ 
tions  (13)  and  (14)  by  -8il/f(kJ,t)l(l+8)  and 
-(/£(/:, /,*)/(  1  +  5),  respectively,  where  the  superscript 
asterisk  represents  a  complex  conjugate,  add  the  two  re¬ 
sults  together,  take  the  real  part,  assume  a  statistically 
steady  state,  and  take  averages  over  time  (so  that  the  dldt 
terms  drop  out),  to  obtain  the  depth-averaged  spectral 
energy  equation 

+  0  +  T’aPE^’O  +  ^forcing(^) 

+  7'fncio„(^0+7'viscosity(/c,/)  =  0.  (15) 

In  (15)  the  spectral  transfers  of  upper  and  lower  layer 
kinetic  energy  are 

TKEA(k,f)  =  Re  /,  l' 0). 


rKa2(M  =  Rc(^pt(kj,t) 

x7(*2,~vV2XM,o).  (17) 

respectively,  where  Re  denotes  the  real  part  of  a  com¬ 
plex  number,  and  the  bracket  operator  ( )  represents  a 
time  average.  The  spectral  transfer  of  available  potential 
energy  is 

TAPE(k,i)  =  Re((1;-;2L^r(M,0 


X/(</q,</f2)(M,0y-  (18) 

The  spectral  transfer  due  to  the  imposed  mean  flow 
forcing  is 

7W 

+  Re(i 


and 


(16) 


(19) 
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the  spectral  transfer  due  to  bottom  Ekman  friction  is 

WM  =  -r^*2*2 (20) 

and  the  spectral  transfer  due  to  eddy  viscosity  is 

T'viscosi.y^/)  =  ~T 

-■^vK*(%{k,l,W2(k,l,t)).  (21) 

As  noted  previously,  the  effects  of  the  wavenumber  fil¬ 
ter  on  the  energy  budgets  will  be  calculated  as  a  residual. 
Because  of  the  filter,  in  practice,  the  right-hand  side  of 
(15)  is  not  identically  zero. 

Spectral  fluxes  arc  defined  as  semi-infinite  integrals  of 
the  spectral  transfers.  For  example 

nKE1(/0=  [f  TKEl(k,l)dkdl,  (22) 

JJki+P^K2 

and  more  generally, 

n,ermW=ffJ  „  Tlenn(k,l)dkdl,  (23) 

J  Jk^+P^K2 

where  the  subscript  “term”  could  denote  “KE,  1,”  “KE,  2,” 
“APE,”  “forcing,”  “friction,”  or  “viscosity.”  Following  the 
terminology  of  Scott  and  Arbic  (2007),  we  will  call  the 
budget  obtained  by  summing  the  spectral  fluxes  II(tf)  [as  in 
Eq.  (15),  but  with  every  “7”  replaced  by  a  “II”]  the  “in¬ 
tegral  energy  budget.”  Note  that  since  we  are  using  the 
depth-averaged  energy  equation  there  are  no  spectral 
transfers  or  fluxes  representing  energy  transfer  between 
layers. 

c.  Impact  of  viscosity  on  idealized  model 

Our  first  look  at  the  spectral  fluxes  I1ke.i(^0  of  upper- 
layer  kinetic  energy  in  the  QG  model  is  displayed  in 
Fig.  3.  Fluxes  are  shown  for  F ^  —  0.4  simulations  with  four 
different  nondimensional  v  values  (0,  0.01,  0.1,  and  1). 
Figure  3a  presents  the  fluxes  normalized  by  ( [u\  -  U2)3ILd, 
while  in  Fig.  3b  all  spectral  flux  curves  are  normalized  by 
their  minimum  (negative)  n(/f)  value  at  larger  scales,  so 
that  the  minimum  values  in  the  normalized  curves  are 
always  -1.  The  spectral  flux  nKEa(/C)  from  the  simula¬ 
tion  with  Fl  ~  0  and  zero  viscosity — with  only  a  wave- 
number  filter  present  to  absorb  the  enstrophy  cascade — is 
represented  by  black  curves  in  Figs.  3a  and  3b.  There  is  a 
vigorous  inverse  cascade  at  scales  near  Ld— both  scales 
slightly  larger  than  Ld  and  scales  slightly  smaller  than 
Ld — but  very  little  forward  cascade  at  small  scales.  The 
small  forward  cascade  at  small  scales  is  consistent  with  the 


(a) 


Fig.  3.  Spectral  kinetic  energy  flux  Hke.K^O  *n  upper  layer  of 
two-layer  QG  model,  for  low-resolution  simulation  with  zero  vis¬ 
cosity,  and  simulations  with  varying  nondimensional  viscosity  v. 
The  low-resolution  experiment  was  run  in  a  larger  domain  with  the 
same  number  of  gridpoints  as  in  the  other  experiments,  such  that 
scales  near  Ld  were  not  as  well  resolved.  In  all  simulations  shown 
above  the  nondimensional  bottom  friction  FL  is  0.4.  In  (a),  the 
fluxes  are  normalized  by  (i7i  -  Ui)3ILd.  In  (b),  all  flux  curves  are 
normalized  so  that  the  minimum  value  in  the  negative  lobe  is  -1. 

fact  that  almost  all  of  the  energy  dissipation  is  accom¬ 
plished  by  bottom  drag  in  the  v  =  0  simulation,  and  very 
little  is  accomplished  by  the  wavenumber  filter. 

In  the  v  =  0.01  experiment  (cyan  curves)  there  is 
a  small  forward  cascade.  In  the  experiments  with  non- 
dimensional  rvalues  of  0.1  and  1  (blue  and  green  curves, 
respectively),  there  is  a  substantial  forward  cascade 
(positive  lobe)  at  smaller  scales,  in  qualitative  agree¬ 
ment  with  the  shape  seen  in  the  observations  (Fig.  I).2 
At  first  glance  this  suggests  that  the  forward  cascades 
seen  in  altimetry  data  could  be  robust  and  due  to  pro¬ 
cesses  that  can  roughly  be  characterized  with  a  hori¬ 
zontal  eddy  viscosity.  We  therefore  ask  whether  0.1  is 
a  plausible  value  for  the  nondimensional  v  in  the  mid¬ 
ocean.  If  we  take  the  dimensional  horizontal  eddy 
viscosity  to  equal  50  m2  s”1,  as  argued  by  Polzin  (2008, 


2  Note  that  the  forward  cascade  is  present  in  both  the  barotropic 

and  baroclinic  modes  (not  shown). 


February  2013 


ARBIC  ET  AL. 


291 


2010),  together  with  typical  midlatitude  open-ocean 
gyre  values  of  Ld  =  50  km  (Richman  et  al.  1977)  and 
u\-U2  =  1  cms"1  (Stammer  1997),  we  obtain  a  non- 
dimensional  v  of  0.1,  suggesting  that  this  may  indeed 
be  a  plausible  value.  (Note,  however,  that  the  non- 
dimensional  v  values  inferred  later  from  the  realistic 
model  are  somewhat  smaller.)  The  spectral  fluxes 
computed  from  experiments  with  larger  nondimen- 
sional  v  values  display  a  shift  of  the  zero  crossings  to 
lower  wavenumbers  (seen  more  easily  in  Fig.  3b).  In 
the  v  -  0.1  experiment  the  zero  crossing  lies  very  close 
to  the  deformation  wavenumber,  as  in  the  observa¬ 
tional  n(AT)  curve  shown  in  Fig.  1.  In  the  v  =  0  experi¬ 
ment  there  is  some  negative  flux  (inverse  cascade) 
occurring  for  wavenumbers  higher  than  the  deformation 
wavenumber. 

The  magenta  curves  in  Figs.  3a  and  3b  denote 
nKE4(/0  computed  from  a  v  -  0,  FL  =  0.4  simulation  in 
which  the  number  of  gridpoints  is  the  same — 2562 — but 
in  which  the  domain  is  4  times  larger  on  each  side.  The 
horizontal  grid  resolution  of  this  run  is  therefore  lower. 
In  the  low-resolution  run  the  spectral  fluxes  point  to 
a  forward  cascade  of  energy  at  small  scales  toward  yet 
smaller  scales,  and  the  subgrid-scale  dissipation  ac¬ 
complishes  a  substantial  fraction  (about  14%)  of  the 
energy  dissipation.  However,  since  this  cascade  is  not 
present  in  the  higher-resolution  v  =  0  simulation,  the 
forward  flux  is  merely  a  result  of  insufficent  horizontal 
resolution.  Note  that  lower  horizontal  resolution  is 
typical  of  many  QG  turbulence  studies,  for  example, 
Larichev  and  Held  (1995),  since  they  often  focus  on  the 
inverse  cascade  and  utilize  large  domains  to  give  the 
inverse  cascade  room  to  proceed.  We  will  return  to  this 
theme,  of  the  effects  of  insufficient  horizontal  resolu¬ 
tion,  shortly.  In  the  remainder  of  the  discussion  in  this 
section  we  return  to  experiments  performed  at  higher 
horizontal  resolution,  though  in  some  analyses  the 
model  results  are  filtered  after  the  run  is  completed. 

Horizontal  eddy  viscosity  significantly  alters  the  QG 
model  energy  budget.  This  can  be  seen  in  Fig.  4,  which 
displays  the  depth-integrated  integral  energy  budget  in 
the  simulation  without  viscosity  (Fig.  4a;  black  curves  in 
Fig.  3)  and  in  the  simulation  with  parameters  otherwise 
equal,  but  with  a  nondimcnsional  v  value  of  0.1  (Fig.  4b; 
blue  curves  in  Fig.  3).  In  the  simulation  without  viscosity, 
bottom  drag  accomplishes  nearly  all  of  the  energy  dis¬ 
sipation.  In  the  simulations  with  v  =  0.1,  about  half  of 
the  energy  dissipation  is  accomplished  by  viscosity.  This 
is  in  reasonable  agreement  with  a  regional  mesoscale 
energy  budget  derived  from  observations  (Polzin  2010), 
which  suggests  that  eddy  viscosity  and  bottom  drag 
dissipate  similar  amounts  of  energy.  Note  that  in  Fig.  4 
the  upper-  and  lower-layer  contributions  to  the  forcing, 


(a)  v=0 


FlG.  4.  Integral  energy  budgets — budgets  of  spectral  fluxes 
U(K) — of  two-layer  QG  simulations  with  FL  =  0.4  and  non- 
dimensional  horizontal  eddy  viscosity  values  v  of  (a)  0  and  (b)  0.1. 
In  both  subplots  the  cyan  line  denotes  the  residual  of  the  other 
terms,  assumed  dominated  by  the  wavenumber  filter.  All  terms 
normalized  by  {u\  -  UiflLj. 

viscosity,  and  kinetic  energy  fluxes  have  been  added 
together;  the  upper  layer  dominates  these  sums. 

Next,  we  briefly  discuss  the  sensitivity  of  QG  turbu¬ 
lence  to  viscosity.  Figure  5  displays  various  eddy  statis¬ 
tics  as  a  function  of  nondimensional  v .  We  use  three 
different  values  of  FL  (0.2, 0.4,  and  0.8).  In  the  first  five 
subplots  of  Fig.  5  extra  horizontal  lines  denote  rough 
ranges  of  values  observed  at  different  locations  in  the 
ocean  (see  also  Table  1;  note  that  the  range  of  observed 
Ftol  values — i.e.  the  ratio  of  peak  forward  to  peak  in¬ 
verse  fluxes-is  taken  from  the  regional  AVISO  calcu¬ 
lations  to  be  described  in  section  4).  Figures  5a-d  display 
relatively  flat  sensitivity  to  v  until  v  takes  on  values  of 
order  one  and  larger.  The  large  v  regime  is  not  likely  to 
be  relevant  to  the  ocean  but  is  included  for  the  sake  of 
completeness  in  the  sensitivity  study.  For  v  of  order  one 
and  smaller,  the  surface  eddy  kinetic  energy  behaves  as 
one  might  expect — it  generally  decreases  with  in¬ 
creasing  v  and  with  increasing  FL  (Fig.  5a).  The  length 
scale  L\  of  upper  layer  kinetic  energy,  defined  by 
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<c>  keWkebt 


(b)  LAd 


(e)  Peak  forward  to  inverse  flux 


10  10  10 
Nondimensional  v 


(f)  Energy  dissipation  by  viscosity/production 


Nondimensional  v 


FlG.  5.  Various  time-averaged  eddy  statistics  in  two-layer  QG  turbulence  model,  vs  nondimensional  v  (*  axes).  The 
jc  axes  are  log  axes;  experiments  with  v  =  0  are  assigned  to  v  *  1CT3  so  that  they  can  be  placed  onto  these  log  axes. 
Three  different  nondimensional  bottom  friction  ( FL )  values  are  used;  see  legend  in  subplot  b.  Extra  horizontal  lines 
denote  rough  ranges  of  values  inferred  from  observations  (see  text  and  Table  1).  (a)  Surface  eddy  kinetic  energy 
(KE)0.5(uf  +  vf),  normalized  by  mean  flow  kinetic  energy  0.5(wi  —  t/2)2(b)  Ratio  of  upper  layer  length  scale  L\  (see 
text)  to  Ld.  (c)  Ratio  KEbcVKEbt  of  baroclinic  to  barotropic  kinetic  energy,  (d)  Ratio  of  time-averaged  squared 
upper-layer  velocities  {u\+  v\)  to  time-averaged  squared  lower  layer  velocities  (u\  +  1^).  (e)  Ratio  Ftol  of  peak 
forward  to  inverse  flux  [see  Eq.  (1)].  (f)  Ratio  Eu  of  energy  dissipation  by  viscosity  to  total  energy  production. 
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Table  1.  Various  time-averaged  eddy  statistics  in  two-layer  QG  turbulence  model,  in  selected  runs  with  different  values  of  v  and  FL  as 
given.  The  v  =  0,  FL  =*  0.4  simulation  is  the  “nominal  solution”  in  Arbic  and  Flierl  (2004),  which  they  argued  compared  well  to  ocean 
observations.  See  text  and  caption  to  Fig.  5  for  definitions  of  symbols  used  in  this  Table.  Rough  ranges  of  observed  values  are  also  given. 
Sources  for  the  ranges  of  observed  values  are  given  by  superscripts  below. 


Surface  KE 

u_ 

Lj 

KEbc 

KEbt 

(“i +ui> 

(«!+”!> 

Ftol 

Ev 

Observations 

10-100“ 

1-3  b 

0.5-3c 

10-100* 

0.17-0.31* 

— 

v  -  0,  Ft  =  0.4 

21 

1.1 

1.5 

31 

0.012 

0 

v  =  0.05,  Fl  =  0.4 

13 

1.4 

1.6 

27 

0.077 

0.32 

p  =  0.1,  Fl  =  0.4 

9.3 

1.5 

1.8 

29 

0.13 

0.45 

v  =  0.05,  Fl  =  0.2 

29 

U 

0.65 

9.9 

0.063 

0.28 

v  =  0.1,  Fl  =  0.2 

20 

1.6 

0.75 

11 

0.10 

0.41 

v  =  0.2,  Fl  =  0.2 

13 

1.7 

0.98 

12 

0.14 

0.57 

v  =  0.01,  Fl  =  0.8 

12 

1.3 

2.8 

77 

0.034 

0.14 

•  Wunsch  (2001). 

b  Stammer  (1997);  see  also  Fig.  9b  of  Arbic  et  al.  (2007)  for  observed  eddy  length  scales  estimated  using  the  same  method  used  here  to 
compute  length  scales  from  the  QG  models. 
c  Wunsch  (1997). 
d  Schmitz  (1996). 
c  Calculations  in  section  4. 


_  -i  -1 

:  +  PE(k,  /)  dk  dl 

-  ,  (24) 

E(k,  [)  dk  dl 


where  Vk2  +  P  =  K  is  isotropic  wavenumber,  and  E(k>  l) 
is  the  wavenumber  spectrum  of  upper  layer  kinetic  en¬ 
ergy,  is  plotted  in  Fig.  5b.  Increasing  nondimensional  v 
leads  to  increasing  eddy  length  scales,  in  keeping  with 
the  smoothing  effect  of  viscosity  seen  in  Fig.  2.  For  large 
v  values  the  flow  becomes  increasingly  baroclinic  (Fig. 
5c)  and  increasingly  surface  intensified  (Fig.  5d).  The 
degree  of  surface  intensification  increases  as  the  bottom 
friction  Ft  increases,  and  is  also  significantly  impacted 
by  the  value  of  8  (Arbic  and  Flierl  2004).  As  non- 
dimensional  values  of  v  increase,  the  Ftol  values  gener¬ 
ally  increase  (Fig.  5e).  Accordingly,  the  ratio  Ev  of  energy 
dissipated  by  viscosity  to  total  energy  production — where 
both  are  taken  from  Eq.  (11) — increases  steadily  with 
increasing  viscosity  v ,  approaching  unity  for  nondimen¬ 
sional  v  values  of  order  one  or  larger.  Figures  5e  and  5f 
demonstrate  that  even  v  values  as  small  as  0.01  yield 
a  nonnegligible  forward  kinetic  energy  cascade  and  vis¬ 
cous  dissipation. 

As  seen  in  Figs.  5a-d,  the  amplitude,  horizontal  scale, 
and  vertical  structure  of  the  eddy  kinetic  energy  in  the 
QG  model  generally  lie  outside  of  observed  ranges  for 
nondimensional  v  values  larger  than  order  one,  and  in¬ 
side  or  at  least  near  the  observed  range  for  v  values  of 
order  one  and  smaller.  As  seen  in  Fig.  5e,  the  Ftol  values 
lie  near  the  range  seen  in  present-generation  altimeter 
data  only  for  v  values  of  about  0.1  or  larger.  However, 
the  observed  range  in  Fig.  5e  may  be  too  large  for  rea¬ 
sons  to  be  discussed  shortly.  Lowering  the  observed  Ftol 


values  would  imply  that  the  small  v  QG  experiments 
match  the  Ftol  values  in  the  data  more  closely. 

To  complement  Fig.  5,  Table  1  displays  values  of  the 
eddy  statistics  for  selected  QG  experiments,  in  which  the 
statistics  in  Figs.  5a-d  lie  close  to,  or  inside,  the  ranges  of 
observed  values.  Taken  together,  Fig.  5  and  Table  1  in¬ 
dicate  that  it  is  possible  for  the  amplitudes,  horizontal 
scales,  and  vertical  structure  of  eddies  in  QG  turbulence 
simulations  to  match  values  in  observations  reasonably 
well,  while  still  having  a  nonnegligible  forward  cascade 
and  viscous  dissipation. 

d.  Effect  of  filtering  on  spectral  fluxes  in  idealized 
model 

Next  we  investigate  the  effect  of  filtering  the  fields 
on  the  computations  of  ITke.i(^)  from  idealized  QG 
model  output.  We  are  particularly  interested  in  simu¬ 
lations  that  do  not  display  a  forward  cascade  in  the  un¬ 
filtered  case,  as  we  wish  to  discern  whether  filtering  will 
yield  an  artificial  forward  cascade.  Note  that  filtering  is 
performed  on  the  model  output  after  the  model  has  been 
run,  not  during  the  run  itself.  Figure  6  shows  ITke,i(^) 
computed  from  the  v  =  0,  FL  —  0.4  QG  simulation 
(black  curves  in  Fig.  3),  and  from  filtered  versions  of  this 
simulation.  The  spatial  filter  goes  to  zero  6  grid  points 
from  the  central  grid  point,  meaning  it  has  a  “full-width 
half  maximum”  (width  at  the  half-power  point)  of  3.4 
grid  points,  or  a  wavelength  of  0.83 L,/.  The  time  filter  is 
a  Blackman  filter  which  goes  to  zero  4  units  of  non- 
dimensional  time  Ldi(u\  -  uf)  from  the  central  time. 
Filtering  in  space,  or  in  time,  diminishes  the  negative 
lobe  of  the  spectral  flux  (inverse  cascade)  substantially. 
Filtering  in  both  space  and  time  further  diminishes  the 
negative  lobe.  Widening  the  time  filter  to  9  units  of 
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FIG.  6.  Spectral  flux  I1ke,i(JQ  of  upper-layer  kinetic  energy  com¬ 
puted  from  the  v  -  0,  FL  =  0.4  two-layer  QG  simulation  (“un¬ 
filtered”)  and  from  filtered  versions  of  this  simulation.  See  text  for 
descriptions  of  filters  used.  All  fluxes  normalized  by  (wi  -  uiflLd- 

nondimensional  time  diminishes  the  inverse  cascade 
even  more.  Most  importantly,  we  see  that  a  simulation 
with  v  —  0  can  see  a  substantial  forward  flux  at  smaller 
scales  if  the  </q  field  is  first  spatially  and/or  temporally 
filtered  before  FIke,i(^0  is  computed,  whereas  little 
forward  flux  is  seen  in  the  fluxes  computed  from  un¬ 
filtered  v  -  0  model  output.  This  suggests  that  spectral 
fluxes  computed  from  the  present-day  AVISO  gridded 
altimeter  product,  which  is  highly  filtered,  could  display 
a  forward  cascade  as  an  artifact  of  insufficient  horizontal 
resolution.  In  contrast,  the  existence  of  the  inverse  cas¬ 
cade  is  robust  to  the  filtering  exercise,  though  the  mag¬ 
nitude  of  the  inverse  cascade  is  not.  Finally,  we  note  that 
filtering,  like  an  increase  of  viscosity,  tends  to  shift  the 
zero  crossings  of  the  nKE,i(/Q  curves  toward  lower 
wavenumbers. 

4.  Results  from  the  realistic  model  and  gridded 
satellite  altimeter  data 

a.  Description  of  realistic  model 

We  use  three  full  years  of  output  (2001-03)  from  the 
forward  (nonassimilative)  version  of  the  Naval  Research 


Laboratory  (NRL)  Layered  Ocean  Model  (NLOM). 
Data-assimilative  versions  of  NLOM  are  in  current  use  as 
a  U.S.  Navy  operational  model.  Shriver  et  al.  (2007)  and 
references  therein  describe  the  NLOM  grid,  wind  forcing, 
and  other  model  details.  The  model  horizontal  resolution 
is  V320  in  latitude  and  45/1024°  in  longitude,  on  a  model 
grid  extending  from  72°S  to  65°N.  For  simplicity,  we  refer 
to  the  simulations  as  having  V320  resolution.  NLOM  is 
based  on  the  primitive  equation  model  of  Hurlburt  and 
Thompson  (1980),  but  with  greatly  expanded  capability 
(Wallcraft  et  al.  2003).  It  has  6  dynamical  layers  and 
a  bulk  mixed  layer. 

Strengths  of  NLOM  for  the  current  study  include  the 
simplicity  of  the  subgrid-scale  closures  it  employs  rela¬ 
tive  to  the  closures  employed  in  other  widely  used  re¬ 
alistic  models.  In  the  NLOM  momentum  equations,  the 
only  subgrid  scale  dissipation  utilized  is  a  horizontal 
eddy  viscosity  of  20  m2  s-1,  comparable  to  the  values 
inferred  from  observations  by  Polzin  (2010).  A  hori¬ 
zontal  eddy  diffusivity  of  equal  strength,  acting  on 
density,  is  also  employed  in  the  model.  In  addition, 
horizontal  biharmonic  diffusivities  of  2.5  X  107  m4  s_1 
are  employed  on  density  and  layer  thickness.  The  only 
vertical  (between  layer)  process  available  in  NLOM  is 
interfacial  friction,  which  is  not  turned  on  in  the  simu¬ 
lations  examined  for  this  paper. 

b.  Satellite  altimeter  data 

We  use  840  snapshots  of  the  two-satellite  AVISO  V30 
Mercator  grid  “reference”  product  (Le  Traon  et  al. 
1998;  Ducet  et  al.  2000),  beginning  with  14  October  1992 
and  ending  with  12  November  2008.  The  AVISO  record 
we  utilize  is  much  longer  than  the  NLOM  record  (16 
versus  3  years)  but  has  a  much  lower  temporal  sampling 
rate  (7  days  versus  6  h).  The  spatial  resolution  is  also 
much  lower  (V30  versus  V320). 

c.  Regions  used  in  analyses 

For  both  the  realistic  model  and  satellite  altimeter 
analyses,  we  compute  spectral  fluxes  over  six  regions.  The 
regions  are  shown  in  Fig.  7,  against  backdrops  of  SSH  and 
SSH  anomaly  snapshots  from  the  model  and  altimeter 
data,  respectively.  Two  of  the  regions — “midlatitude 
southeast  Pacific”  and  “high  latitude  southeast  Pacific” — 
are  relatively  quiescent.  Four  of  the  regions — “Agulhas,” 
off  the  tip  of  southern  Africa,  “Malvinas,”  in  the  western 
South  Atlantic,  “Gulf  Stream,”  in  the  western  North 
Atlantic,  and  “Kuroshio,”  in  the  western  North  Pacific — 
contain  strong  western  boundary  currents  and  mesoscale 
eddy  activity,  as  seen  in  Fig.  7.  The  correspondence  be¬ 
tween  regions  in  the  AVISO  versus  NLOM  analysis  is  not 
exact,  because  of  the  lack  of  grid  points  in  shallow  waters 
in  NLOM  (see  Shriver  et  al.  2007  and  references  therein). 
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FiG.  7.  Regions  used  to  compute  spectral  fluxes  riKE,i(^0  of  surface  ocean  geostrophic  ki¬ 
netic  energy  in  (a)  Vyf  NLOM  output  and  (b)  V3°  Mercator  AVISO  gridded  satellite  altimeter 
data.  The  NLOM  regions  are  highlighted  against  the  15  Aug  2002  snapshot  of  sea  surface  height 
(cm)  in  the  model,  while  the  AVISO  regions  are  highlighted  against  the  20  Dec  2006  snapshot 
of  sea  surface  height  anomaly  (cm)  in  AVISO.  For  the  sake  of  comparison,  for  this  figure  the 
AVISO  data  was  interpolated  onto  a  V?  regular  latitude-longitude  grid  having  the  same  lat¬ 
itudinal  range  as  the  NLOM  grid. 


The  misfit  between  the  boundaries  of  the  NLOM  versus 
AVISO  regions  is  largest  in  the  Agulhas. 

d.  Spectral  fluxes  in  realistic  model  and  gridded 
altimeter  data 

We  now  describe  the  spectral  flux  calculations  per¬ 
formed  with  output  from  the  realistic  model  (NLOM) 
and  from  gridded  satellite  altimeter  data  (AVISO).  The 
NLOM  sea  surface  heights  are  first  low-pass  filtered  tn 
time  to  remove  motions  with  periods  of  3  days  or  less. 
This  is  done  because  high-frequency  motions  cannot  be 
geostrophic.  For  both  the  NLOM  and  AVISO  analyses, 
we  compute  \p\  =  grj//,  where  17  is  the  perturbation  SSH, 
g  —  9.8  m  s~2  is  gravitational  acceleration,  and /is  the 
Coriolis  parameter,  which  varies  as  the  sine  of  the 
gridpoint  latitude  (Vallis  2006).  Then  for  each  of  the  six 
regions  shown  in  Fig.  7,  each  snapshot  of  ipi  is  detrended 
in  space,  with  a  two-dimensional  least  squares  fit.  Fol¬ 
lowing  this,  the  ip  1  field  is  tapered  in  space  with  a  two- 
dimensional  function  constructed  of  nine  overlapping 
Hanning  windows  in  each  spatial  direction,  where  each 
window  has  a  width  one-fifth  of  the  domain  extent  in 
that  direction.  This  type  of  detrending  and  tapering  is 


standard  in  computations  of  Fourier  transforms  of  data 
in  nonperiodic  domains  (Priestley  1981).  The  sensitivity 
to  tapering  will  be  discussed  shortly.  Fourier  transforms 
arc  repeatedly  utilized  to  compute  derivatives  and 
hence  the  term  /(i/q,  V2^).  For  simplicity  we  base  the 
increment  lengths  of  the  wavenumbers  in  the  east-west 
direction  on  the  length  of  the  increments  at  the  central 
latitude  of  the  box.  In  other  words,  we  do  not  account 
for  the  convergence  of  meridians  in  our  spectral 
calculations.  As  in  the  QG  calculations,  the  transfer 
T — the  real  part  of  where  the  asterisk 

indicates  the  complex  conjugate — is  computed  for  all 
snapshots,  and  then  averaged  over  the  snapshots.  An 
integration  over  the  resulting  time-averaged  T  values  for 
all  wavenumbers  greater  than  a  given  wavenumber  K 
yields  the  spectral  flux  IIke,i(^0  for  that  wavenumber. 

Figures  8-10  display  the  spectral  flux  of  surface  ocean 
geostrophic  kinetic  energy  nKEtl(AT)  computed  from  the 
NLOM  simulation,  alongside  fluxes  computed  from 
AVISO  gridded  altimeter  data,  for  the  six  regions  of 
interest.  Note  that  in  these  more  realistic  computations 
nKE.i(/0  does  not  include  the  factor  of  5/(1  4-  5)  uti¬ 
lized  in  section  3.  This  factor  is  dropped  here  for 
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(a)  Midiatftude  SE  Pacific 


Fig.  8.  Spectral  flux  nKE,i(^0  °f  surface  ocean  geostrophic 
kinetic  energy  in  (a)  midlatitude  southeast  Pacific  and  (b)  high- 
latitude  southeast  Pacific,  computed  from  AVISO  (cyan  curves) 
and  from  four  versions  of  realistic  NLOM  simulation  output:  un¬ 
filtered  (black  curves),  filtered  in  space  (red  curves),  filtered  in  time 
(blue  curves),  and  filtered  in  space  and  time  (green  curves).  See 
text  for  details.  Vertical  dashed  lines  denote  deformation  wave- 
number  1  ILd,  where  the  deformation  radius  Ld  is  taken  from 
Chelton  et  al.  (1998). 

simplicity,  since  oceanic  stratification  varies  from  one 
location  to  another. 

The  black  curves  in  Figs.  8-10  show  the  “unfiltered” 
NLOM  fluxes — that  is,  the  fluxes  computed  from 
NLOM  output  that  is  unadalterated  except  for  the  ma¬ 
nipulations  described  in  the  first  paragraph  of  this  sub¬ 
section.  Consistent  with  the  results  of  Schlosser  and 
Eden  (2007)  and  Tulloch  et  al.  (2011),  the  spectral  flux 
computed  from  a  realistic  high-resolution  model  has  a 
prominent  negative  lobe,  indicating  an  inverse  cascade, 
at  larger  scales.  The  deformation  wavenumbers  1  !Ld, 
where  Ld  is  the  first  baroclinic  mode  deformation  radius, 
are  shown  as  vertical  dashed  lines  on  Figs.  8-10.  The  L d 
values  are  taken  from  the  atlas  of  Chelton  et  al.  (1998). 
For  each  of  the  six  regions,  we  choose  the  Ld  value  at  the 
atlas  gridpoint  lying  closest  to  the  center  of  the  AVISO 
region  in  question.  In  all  six  regions,  NLOM  displays 
a  significant  inverse  cascade  (negative  spectral  fluxes) 
over  a  band  of  wavenumbers  higher  than  the  deforma¬ 
tion  wavenumber,  consistent  with  the  behavior  seen  in 


(a)  Agulhas 


(b)  Malvinas 


Fig.  9.  As  in  Fig.  8,  but  for  (a)  Agulhas  and  (b)  Malvinas. 

the  QG  results  with  zero  or  low  viscosity  (i.e.,  black  and 
cyan  curves  tn  Fig.  3). 

Table  2  lists  the  estimated  values  of  nondimensional  v , 
and  of  Ftol  (the  ratio  of  peak  forward  to  peak  inverse 
fluxes),  for  the  six  regions  of  NLOM  output.  Recall  that 
nondimensional  v  in  the  two-layer  QG  model  equals  the 
dimensional  viscosity  divided  by  the  product  of  Ld  and 
the  mean  shear  (u\  -  Hi).  The  NLOM  dimensional  v 
value  was  given  in  section  4a.  For  each  of  the  six  re¬ 
gions,  the  deformation  radius  Ld  is  obtained  from  the 
Chelton  et  al.  (1998)  dataset  as  described  above.  To 
estimate  the  NLOM  mean  shear  (the  NLOM  equiva¬ 
lent  of  U\  —  U2  in  the  two-layer  QG  model),  we  com¬ 
puted  the  time-averaged  velocity  components  Hi ,  1/6,  vi , 
and  V6  at  all  model  gridpoints  from  year  2003  of  the 
NLOM  simulation,  where  here  subscripts  1  and  6  re¬ 
spectively  denote  the  top  and  bottom  layers  of  the  six 
layers  in  NLOM.3  Then  over  each  of  our  six  regions  of 
interest,  we  computed  the  area-averaged  RMS  shear  as 


3  Top-  and  bottom-layer  NLOM  velocities  were  not  saved  over 
all  three  of  the  years  2001-03  analyzed  in  this  paper.  Note  that  we 
are  analyzing  NLOM  runs  performed  some  time  ago,  for  purposes 
other  than  those  envisioned  for  the  present  paper. 
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(a)  Guff  Stream 


FlG.  10.  As  in  Fig.  8,  but  for  (a)  Gulf  Stream  and  (b)  Kuroshio. 


where  dA  represents  an  element  of  area.  Table  2  dem¬ 
onstrates  that  the  NLOM  nondimensional  v  values  are 
of  order  0.01,  suggesting  that  our  earlier  estimates  of  v  = 
0.1  for  the  ocean  may  have  been  too  large.  Table  2  shows 
that  the  Ftol  values  generally  decrease  with  decreasing 
values  of  nondimensional  v ,  in  qualitative  agreement 
with  the  QG  results.  To  make  a  more  quantitative 
comparison,  we  note  that  the  v  —  0.01  QG  experiment 
listed  in  Table  1  has  an  Ftol  value  of  0.034,  which  is  of 
the  same  order  as  the  Ftol  values  in  the  first  four  regions 
listed  in  Table  2,  all  of  which  have  nondimensional  v 
values  of  the  same  order  as  0.01. 

Figures  8-10  also  display  the  spectral  fluxes  computed 
from  NLOM  output  that  has  been  filtered  in  space,  time, 
or  both  space  and  time,  in  an  attempt  to  mimic  the  fil¬ 
tering  inherent  in  the  construction  of  AVISO  gridded 
altimeter  data.  As  in  the  results  computed  from  filtered 
output  of  the  idealized  model,  the  filters  act  on  model 
output,  not  on  the  model  itself  as  it  is  running.  We 
choose  to  use  Blackman  filters  in  both  space  and  time. 
Since  the  AVISO  product  is  put  out  every  seven  days, 


Table  2.  Computed  nondimensional  viscosities  v  and  Ftol 
values  (see  text)  for  six  regions  of  NLOM  output.  Results  listed  in 
order  of  descending  v. 


Region 

Nondimensional  v 

Ftol 

Midlatitude  southeast  Pacific 

0.015 

0.050 

High-latitude  southeast  Pacific 

0.011 

0.024 

Malvinas 

0.0076 

0.0076 

Kuroshio 

0.0039 

0.014 

Gulf  Stream 

0.0038 

0.0071 

Agulhas 

0.0029 

0.0050 

on  a  V30  grid,  we  choose  a  temporal  Blackman  filter 
which  goes  to  zero  seven  days  from  the  central  time,  and 
a  spatial  Blackman  filter  which  goes  to  zero  at  a  distance 
equal  to  37  km  (V30  of  latitude)  from  the  central  grid- 
point.  Thus  the  chosen  Blackman  filters  have  length  and 
time  scales  comparable  to  those  used  in  constructing 
the  AVISO  product,  although  the  analogy  is  of  course 
rough.  The  effect  on  the  computed  spectral  fluxes  is 
similar  to  the  effect  seen  in  the  fluxes  computed  from 
filtered  versions  of  the  QG  model  output.  Filtering  in 
space  again  increases  the  estimated  forward  cascade, 
even  in  cases  where  there  was  little  forward  cascade  in 
the  unfiltered  realistic  model  output.  Again  we  see  the 
possibility  that  the  forward  cascade  shown  in  Scott  and 
Wang  (2005)  may  be  an  artifact  of  the  relatively  low 
resolution  of  present-day  gridded  altimeter  data.  Fil¬ 
tering  in  either  space  or  time  diminishes  the  strength  of 
the  inverse  cascade,  and  filtering  in  both  further  di¬ 
minishes  the  inverse  cascade.  Filtering  in  either  space  or 
time  also  tends  to  shift  the  zero-crossings  of  the  spectral 
flux  curves  toward  lower  wavenumbers. 

Figures  8-10  include  the  AVISO  spectral  fluxes, 
which  share  some  characteristics  of  the  filtered  NLOM 
results.  In  general  the  AVISO  forward  fluxes  are  sub¬ 
stantially  larger  than  the  forward  fluxes  in  their  un¬ 
filtered  NLOM  counterparts.  The  AVISO  Ftol  values 
range  from  0.17  to  0.31  over  the  six  regions.  The  AVISO 
fluxes  are  shifted  toward  lower  wavenumbers  relative  to 
the  fluxes  computed  from  unfiltered  NLOM  output. 
Finally,  in  contrast  to  the  unfiltered  NLOM  results,  the 
zero  crossings  of  the  AVISO  fluxes  lie  at  wavenumbers 
equal  to  or  lower  than  the  deformation  wavenumber. 

We  compare  the  strength  of  the  inverse  cascade  in 
AVISO  versus  NLOM  by  computing  the  ratio  of  the 
minimum  of  the  negative  lobe  at  larger  scales  in  the 
AVISO  riKE,i(JQ  curve  to  the  minimum  negative 
nKE,i(A0  value  at  larger  scales  computed  from  unfiltered 
NLOM  output.  This  ratio  is  0.14  in  the  high-latitude 
Southeast  Pacific,  1.1  in  the  Kuroshio,  2.1  in  the  Malvi¬ 
nas,  and  between  0.4-0.6  in  the  other  three  regions.  This 
order  of  magnitude  agreement  in  five  of  the  six  regions  is 
encouraging  considering  that  the  minimum  value  of  the 
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negative  lobes  changes  by  two  orders  of  magnitude 
across  the  six  regions.  The  order  of  magnitude  agree¬ 
ment  bolsters  the  claims  of  Scott  and  Wang  (2005)  that 
the  inverse  cascade  at  larger  scales  is  a  robust  feature  of 
the  oceanic  general  circulation. 

The  computed  spectral  fluxes  are  sensitive  to  the 
choice  of  spatial  windows.  The  subject  of  windows  for 
spectral  computations  in  nonperiodic  data  has  received 
a  great  deal  of  scrutiny  over  the  years  (e.g.,  Priestley 
1981).  The  window  utilized  here,  constucted  from  9 
overlapping  Hanning  windows  in  each  horizontal  di¬ 
rection,  is  close  to  unity  over  most  of  the  domain  and 
drops  smoothly  to  zero  at  the  edges.  We  also  experimented 
with  a  Tukey  window  having  a  ratio  of  taper  to  constant 
sections  equal  to  0.2.  The  spectral  fluxes  computed  from 
NLOM  with  the  Tukey  window  (not  shown  here)  are  very 
nearly  equal  to  those  computed  from  the  9  overlapping 
Hanning  windows  employed  in  this  paper.  Spectral 
fluxes  computed  from  coarser-resolution  AVISO  output 
arc  more  sensitive  to  the  usage  of  Tukey  versus  over¬ 
lapping  Hanning  windows  than  are  fluxes  computed  from 
NLOM  output;  AVISO  fluxes  can  differ  by  up  to  a  factor 
of  about  2  between  the  two  cases,  for  some  wavenumbers. 
However,  even  in  the  AVISO  fluxes  the  wavenumbers  at 
which  forward  cascades  are  seen,  and  at  which  inverse 
cascades  arc  seen,  are  virtually  identical.  Therefore  none 
of  our  qualitative  conclusions  about  the  AVISO  spectral 
fluxes,  or  our  quantitative  conclusions  about  the  NLOM 
spectral  fluxes,  appear  to  be  overly  sensitive  to  the  choice 
of  spatial  window,  as  long  as  the  window  lies  near  one  over 
most  of  the  domain. 

5.  Summary  and  discussion 

In  this  paper  we  have  examined  the  effects  of  hori¬ 
zontal  eddy  viscosity  and  of  horizontal  grid  resolution  on 
geostrophic  turbulence,  with  a  particular  focus  on  the 
spectral  fluxes  FIke.i(K)  °f  surface  ocean  geostrophic 
kinetic  energy  in  the  isotropic  wavenumber  (K)  domain. 
We  have  utilized  outputs  from  idealized  two-layer 
quasigeostrophic  (QG)  turbulence  model  simulations 
and  realistic  global  numerical  model  (NLOM)  simula¬ 
tions  in  this  endeavor,  and  compared  the  model  spectral 
fluxes  to  fluxes  computed  from  present-day  gridded 
satellite  altimeter  products  (AVISO).  We  are  motivated 
by  considerations  of  ocean  energetics  and  by  the  wide¬ 
spread  use  of  horizontal  eddy  viscosity  in  models.  The 
addition  of  viscosity  to  the  QG  model  yields  a  forward 
energy  cascade  at  length  scales  near  the  deformation  ra¬ 
dius  toward  smaller  scales,  as  is  apparently  seen  in  com¬ 
putations  of  ocean  surface  spectral  kinetic  energy  fluxes 
computed  from  AVISO  data.  Viscosity  arguably  repre¬ 
sents  the  physics  of  mesoscale  coupling  to  submesoscale 


motions  and  internal  waves.  Viscosity  can  dissipate  a 
substantial  amount  of  energy  in  the  QG  model.  A  sensi¬ 
tivity  study  in  section  3  of  this  paper  indicates  that  vis¬ 
cosity  impacts  the  statistics  of  eddies  in  QG  turbulence 
models,  but  that  the  model  statistics  can  lie  near  observed 
ranges  as  long  as  the  viscosity  is  not  too  large. 

Insufficient  horizontal  resolution — be  it  the  resolution 
chosen  to  run  the  QG  simulations  at,  or  a  deliberate  fil¬ 
tering  in  space  and/or  time  to  downgrade  the  resolution 
of  the  QG  or  realistic  ocean  model  output  after  the 
simulation  is  completed — can  also  yield  a  forward  cas¬ 
cade  at  smaller  scales.  This  filtering  mimics  the  filtering 
of  along-track  satellite  altimeter  data  onto  the  AVISO 
grids.  It  is  therefore  difficult  to  tell  whether  the  forward 
spectral  kinetic  energy  fluxes  observed  for  small  scales  in 
the  Scott  and  Wang  (2005)  altimeter-based  calculations 
of  nKE,iW  are  physical  or  an  artifact  of  the  smoothing 
inherent  in  constructing  the  AVISO  gridded  altimeter 
product. 

Similar  conclusions  were  reached  by  Boer  and 
Shepherd  (1983)  in  their  pioneering  computation  of 
spectral  kinetic  energy  fluxes  in  atmospheric  data. 
T.  Shepherd  2010  (personal  communication)  notes 
that  “explicit  calculation  of  nonlinear  interactions  is 
inherently  limited  by  the  resolution  of  the  data  set,  espe¬ 
cially  close  to  the  resolution  limit.”  Boer  and  Shepherd’s 
(1983)  attempt  to  account  for  missing  nonlinear  interac¬ 
tions  at  the  limiting  resolutions  of  their  dataset  sub¬ 
stantially  reduced  the  estimated  forward  energy  flux  at 
their  highest  wavenumbers  (see  Fig.  13a  of  their  paper).  In 
addition,  Shepherd  (1987)  found  that  the  crossover  from 
negative  to  positive  flux  occurs  at  different  wavenumbers 
for  different  truncations  of  atmospheric  datasets,  such 
that  the  forward  flux  is  not  a  reliable  feature.  A  simple 
interpretation  of  the  limiting  effect  of  spatial  resolution 
was  given  in  the  2010  personal  communication  with 
Shepherd:  “Following  Fjortoft  (1953),  when  computing 
triads  over  a  finite  spectral  range,  energy  must  cascade 
both  upscale  and  downscale  and  so  an  apparent  downscale 
cascade  is  inevitable  in  such  a  calculation  even  if  there  is 
no  such  cascade  in  nature.  It’s  also  obvious  that  the 
enstrophy  flux  cannot  really  go  to  zero  at  the  truncation 
wavenumber;  this  too  points  to  the  fact  that  resolved 
fluxes  are  strongly  affected  by  truncation.” 

The  proposed  wide-swath  satellite  altimeter  (Fu  and 
Ferrari  2008;  Fu  et  al.  2012)  will  map  sea  surface  heights 
at  much  higher  horizontal  resolution  than  is  currently 
possible.  Energy  cascades  to  small  scales  ultimately  lead 
to  mixing.  Mixing  in  the  upper  ocean  has  a  greater  im¬ 
pact  on  the  oceanic  meridional  overturning  circulation 
than  mixing  in  the  abyssal  ocean  (Scott  and  Marotzke 
2002).  This  adds  importance  to  the  quest  to  understand 
whether  the  forward  fluxes  seen  in  present-generation 
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altimeter  data  represent  real  physics  or  sampling  arti¬ 
facts.  Sea  surface  height  measurements  can  only  be  used 
to  constrain  geostrophic  velocities.  Therefore,  even  data 
from  a  wide-swath  altimeter  can  only  be  used  to  com¬ 
pute  spectral  fluxes  arising  from  the  geostrophic  part  of 
the  flow.  However,  wide-swath  altimeter  measure¬ 
ments  of  sea  surface  height  at  small  scales  will  be  able 
to  constrain  the  height  fields  of  models  at  small  scales, 
and  from  such  models  (e.g.,  Capet  et  al.  2008)  the 
spectral  fluxes  arising  from  ageostrophic  flows  can  be 
computed. 

We  have  shown  that  the  maximum  strength  of  the 
inverse  kinetic  energy  cascade  [i.e.  the  minimum  of  the 
negative  lobes  in  the  spectral  fluxes  nKE1(/Q]  computed 
from  realistic  models  and  altimetry  generally  agree  to 
within  a  factor  of  about  2,  even  as  the  spectral  flux  itself 
varies  by  two  orders  of  magnitude  over  the  six  oceanic 
regions  studied  here.  As  far  as  we  know,  this  is  the  first 
direct  comparison  of  inverse  cascades  computed  in 
realistic  ocean  circulation  models  versus  altimeter 
data.  In  both  the  QG  and  realistic  models,  the  existence 
of  the  negative  spectral  fluxes  (inverse  cascades)  at 
larger  scales  seen  in  current-generation  altimeter  data 
is  robust  to  the  filtering  exercise,  but  the  exact  strength 
of  the  inverse  cascades  is  not  robust  to  filtering.  The 
planned  wide-swath  altimeter  should  therefore  aid  in 
refining  estimates  of  the  strength  of  the  oceanic  inverse 
cascade  as  well  as  constraining  dynamics  at  smaller 
scales. 
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